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Abstract 

The aim of this paper is to apply techniques of symmetry group anal- 
ysis in solving two systems of conservation laws: a model of two strictly 
hyperbolic conservation laws and a zero pressure gas dynamics model, 
which both have no global solution, but whose solution consists of sin- 
gular shock waves. We show that these shock waves are solutions in the 
sense of 1-strong association. Also, we compute all protectable symme- 
try groups and show that they are 1-strongly associated, hence transform 
existing solutions in the sense of 1-strong association into other solutions. 
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The concept of classical symmetry groups offers a large number of possibilities in 
studying differential equations, in particular in constructing explicit solutions 
to linear and nonlinear differential equations or determining and classifying 
invariance properties [161 I17j . In various problems of mathematical physics 
the classical theory turns out to be insufficient, due to singular objects (like 
distributions or discontinuous nonlinearities) which can occur in the equation 
or equations with solutions in a weak sense, i.e. weak solutions (distributional, 
generalized or in the sense of association). Therefore, the methods of classical 
symmetry group analysis of differential equations have been extended to linear 
equations in the class of distributions l2j , as well as to equations involving 
generalized functions [5l [6l [TOl (TTJ [9] . 

The aim of this paper is to apply techniques of symmetry group analysis in 
solving two systems of differential equations given in the form of conservation 
laws. The paper is divided into two parts. Section 1 provides a brief overview 
of the basic definitions and theorems which are going to be used for studying 
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conservation laws. We start by recalling some facts on symmetry group analysis, 
which are in detail carried out in [16] (see also [17]). Then we turn to symme- 
tries in the generalized setting, precisely to associated ones. As we will see later, 
the reason for this lies in the fact that the solutions of the conservation laws we 
consider are, under certain assumptions, shock waves. Lack of space prevents us 
from also giving a short introduction to generalized functions. Therefore, for the 
notations and properties of the Colombeau algebra of generalized functions we 
recommend [Sill] or [T3]; in particular, definitions of generalized step- functions, 
splitted delta functions and m— and m'— singular delta functions are provided 
in [13]. We close the introductory part by a short overview of conservation laws. 
Based on [31 S] [T] [IH] we fix notations and present the general solution of 
the Riemann problem. Motivated by [8] and [13] we proceed in section 2 by 
investigating two systems of conservation laws: a model of two strictly hyper- 
bolic conservation laws which is genuinely nonlinear but for which the Riemann 
problem has no global solution and a zero pressure gas dynamics model which 
is linearly degenerative but for which the Riemann problem also does not have 
global solutions. In both cases singular solutions appear, called singular shock 
waves. We prove that these solutions are solutions in the sense of 1-strong as- 
sociation. After computing all projectable symmetry groups of these systems 
we show that they are 1-strongly associated, hence transform existing solutions 
(given in [S] and [13]) into other solutions. 

1 Introduction 

1.1 Symmetry Groups of Differential Equations 

Let 5 be a system of differential equations: 

A^(x,u) = 0, l<iy<L 

Denote hy X = Mp and U — R"^ the spaces of independent and dependent vari- 
ables with coordinates x = (xi, X2, . • . , Xp) and u = (u^, u^, . . . , u*) respectively. 
Also, denote by M an open subset oi X x U. Identify a function u = f{x) with 
its graph 

Tf = {(x,/(x)) -.xenjcXxU, 

where 51 C X is the domain of /. Let G be a local group of transformations 
acting on M. The transform of F/ by g G G is defined by 

5 ■ T/ = {(i, w) = g • {x, u) : (x, u) € L/}. 
In local coordinates this action is given by 

5 • / = (*s ° {idx X /)) o {Eg o {idx X /))"\ 

where Eg and '^g are smooth function on M , and idx is the identity mapping 
on X. Supposing that Eg does not depend on the dependent variables we get 
projectable action of g on /, i.e. 

g-f={<^go{tdx-^f))oE-g\ (1) 
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1.1 Definition. The symmetry group of the system S is a local transformation 
group G acting on the space of independent and dependent variables with the 
property that whenever u = f{x) is a solution of the system and g ■ f is defined, 
g G G, then u = g ■ f is also a solution of the system. 

The n-th prolonged or n-jet space X x is a space which represents aU 

independent variables, dependent variables and all different partial derivatives 
of dependent variables up to the order n. For the construction of the n-th 
prolonged space we refer to [TB]. Write M*^"^ for a subset of n-jet space 
An arbitrary point in {/("^ will be denoted by w^"^ and its components by u", 
where I < a < q, while J runs over the set of all unordered multi-indices 
J = {ji, ■ ■ -Jk), l<jk<P,0<k<n. 

The n-th prolongation of a function f : X U, denoted by pr^"^^ f ^ is a function 
from X to J/^"-*, which maps x into {d,jf°'{x))a.j, 1 < a < < | J| < n. 
The n-th prolongation of a group G which acts on M C X xU, pr^")G, is again 
a local group of transformations which acts on Af'"^ such that it transforms the 
derivatives of a smooth function u = f{x) into the corresponding derivatives of 
the transformed function u = f{x). For the precise definition see [16] . 
The n-th prolongation of a vector field v on M, pr'^"W, is a vector field on the 
n-jet space M^") with the following property: 

pr(")(exp(?7v))(x,M(")), 

T7 = 

where exp(?7v) is the corresponding local one-parameter group generated by v. 
If 

p a ' d 

i—l a—1 

then we calculate the n-th prolongation of v using the formula: 

pr(")v = v + ^^</);^(x,«("))^, (2) 

where the coefficients 0^(2;, u*^")) are given by 

<i>iix,u^-^) = DjU^-j2eur\ +Ef<.^, (3) 

V i=l / i=l 

uf = du^/dx"^, u°'j^ = du'^/dx^ and Dj denotes a total differential. 

Then the infinitesimal criterion for a system of differential equations reads: 

1.2 Theorem. Let 

A,(a;,u(")) =0, v^l,...,l (4) 

he a system of differential equations of a maximal rank ( meaning that the cor- 
responding Jacohian matrix JA(a^,u'"-') = {dlS.y/dx^,dlS.y/du") is of rank I on 
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the set of all solutions of S, Sa)- If G is a local transformation group acting on 
M CX xU and 

pr(")v(Ai,(x,M("))) = 0, 1^=1,..., I, whenever A(x,u(")) = 0, (5) 

for every infinitesimal generator v of G then G is a symmetry group of 

The condition ^ from this theorem will also be necessary if we additionaly 
suppose that the system (|H) is locally solvable, i.e. at each point (xq, Uq"') € 
there exists a smooth solution u = f{x) of the system, defined in a neighborhood 
of xq, which has the prescribed "initial conditions" Wq"-* — pA'^'' f (xq) ■ (We say 
that a system of differential equations is nondegenerate if at every point of the 
solution set it is both locally solvable and of maximal rank.) 
For later use, we mention here a result which is a consequence of the maximal 
rank condition (imposed on the system ^ in the above theorem) . Namely, un- 
der the conditions of Theorem ll.21 the infinitesimal criterion ([5]) can be replaced 
by the equivalent condition 

I 

pr(")v(A,(a;,ii(")))-^g,^A^(a;,ii(")), z. = l,...,/, (6) 
for functions Qi/^, fi, v = I, . . . ,1 to he determined. 

We finish this short introduction into symmetry groups of differential equations 
by a description of a procedure for calculating symmetry groups of a given 
system S. The procedure consists of the following steps: 

(1) Write the vector field (i.e. infinitesimal generator) in the most general 
form: 

p d d 

i—l a—1 

where and 0^ are functions which should be calculated. 

(2) According to ^ and ^ calculate the corresponding prolongation of v. 

(3) Then apply the infinitesimal criterion (O and equate pr(")v(Aj,(a:, u^"^)) 
with zero. Since those equations must hold on 5'a, eliminate the depen- 
dence of derivatives of u by the equations from the system. After that we 
have the equations which have to be satisfied with respect to x, u and the 
remaining partial derivatives of u. 

(4) After solving these equations we obtain a certain number of partial differ- 
ential equations for and 4>a- 

(5) Compute the and (pa from them, thereby computing vector fields v 
which generate a Lie algebra of infinitesimal symmetries. 

(6) At the end find the corresponding one-parameter symmetry groups as the 
flows of the infinitesimal generators calculated in the previous step. 
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1.2 Symmetry Groups of Weak Solutions 

Next, we look for the symmetries which transform weak solutions of the system 
of PDEs 

A,(x,u(")) = 0, l<iy<l, (7) 

into other weak solutions, mainly associated solutions to ([7]) into other associ- 
ated solutions to d?!) (hence the system ([7]) should be replaced by A^(x,u^"^) « 
0). Such a symmetry group is called symmetry group in the sense of association 
or associated symmetry group for short. The symmetry groups we are interested 
in are projectable in order to avoid the problem of inverting Colombeau func- 
tions. Thus unless explicitly stated otherwise, all symmetry groups are assumed 
to be projectable. Beside this, we need some more assumptions on G (cf. [5]). 
We suppose that a local transformation group G is slowly increasing, uniformly 
for X in compact sets, and analogously for the mapping m*^") i-^ A(a;, u*-"-*). 
Look at the system 

A,(x,u("))«0, 1/ = !,...,/. (8) 
We recall the following definitions from [5]: 

1.3 Definition, u = G Gi^)"^ is a solution of {3) and also as- 
sociated solution to ^ if u = {u^ , ...,u'') has a representative [(uj, uD^] S 
£71/ (f2)'? such that for each test function ip G ^(Q) 

J A^{x,pr'-"'^u,{x))(p{x)dx as l<v<l. (9) 

The set of all associated solutions to (O is denoted by Aa and moreover, the 
set of all u G (^00)' which satisfies ^ with ABa- The symmetry group G of ([7]) 
is called ^-symmetry group if for every u £ A and every G G it follows that 
gr/U G A, whenever g^iU is defined. Beside solution in the sense of association 
we can also define a solution in the sense of strong association. 

1.4 Definition. Let fc G Nq. u= ...,u'^) G Aa{^) (resp. u G ABa{^)) is 

k 

called k-strongly associated or ^.-associated solution to the system ^ if there 
exists a representative [{u\, ...^u1)e\ G SMi^)"^ such that for each B C C'^{Q) 
which is hounded in (11) we have 

lim sup I / A^{x, pr(")Me(a;))(^(a;) dx\ ^ 0, l<i^<l. 

The space of all fc-strongly associated solutions to ([7]) is denoted by AS'}^. Also, 
ABS'}^ := .45^ n Goo- The main role in the calculation of associated, resp. 
/c-strongly associated symmetry groups is played by the theorem which is based 
on the factorization property of the system, derived in In matrix form this 
property is given by 



5 



A(5„(a;, pr^") (5,u)(S,(x, u(x)))) = g(77, pr(")zi(x))A(x, pr(")«(a;)), 

(10) 

where Q : W ^ M'' and W is an open subset of (— J/c'yo) x -W") with {0} x 
M(") C W. 

1.5 Theorem. Let G be a slowly increasing symmetry group of the system 
which admit a global factorization of the form !ilO\}. Then 

(i) if Q depends only on rj, x and u then G is also an ABS symmetry group 

(a) if Q depends only on rj and x then G is also an AS\-symmetry group of 
for each k > 0. Moreover, G is in this case an associated symmetry 
group of ^ as well. 

1.3 Systems of Conservation Laws 

We look at a system of conservation laws in one space dimension: 

< + = 0, 

or written in a shorter (matrix) form: 

ut + {f{u))^ = 0, (11) 

where t > 0, x E M., u = (w^,...,u") is the conserved density and f{u) = 
(/i, . . . , /„) is the flux. Differentiating we obtain a quasilinear system 

Ut + A{u)u^ = 0, (12) 

where A{u) = Df{u) is the Jacobian matrix of/. The systems (fTTj) and (fT2|) are 
equivalent for all smooth solutions u. Otherwise, if u has a jump, the left hand 
side of ()12p contains a product of a discontinuous function with a distributional 
derivative, while (jlip is still well defined in the distributional sense. 
The eigenvalues of the matrix A(u) determine the system of conservation laws 
in the following way: 

1.6 Definition. The system of conservation laws is hyperbolic, resp. strictly 
hyperbolic, if all eigenvalues of the matrix A{u) are real, resp. real and different. 

Suppose that the system (fTTj) is strictly hyperbolic and denote by Xi{u), . . . , 
Xn{u) the eigenvalues of A{u) with Ai(m) < ••• < A„(m). Next, denote by 
li, . . . ,ln and ri, . . . , r„ the corresponding left and right eigenvectors. The eigen- 
value Xi of A is also called the i-th characteristic speed and the pair (A^, r^) the 
i-ih characteristic field of (fTTj) . 
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1.7 Definition. The i-th characteristic field of the system ill]) is called gen- 
uinely nonlinear if 

DX,{u)-r,{u)^0, Vm. 
The i-th characteristic field is called linearly degenerate if 

D\i{u) ■ ri{u) — 0, Vw. 

If the solution of pT|) is a piecewise smooth function u = u{t, x) having a 
discontinuity across a line x — 7(t) with = lim2.^^(t)± u{t, x), then it satisfies 
([T^ outside the 7, while along the line of discontinuities the Rankine-Hugoniot 
conditions holds: 

iu+-u-)j = f{u+)-fiu-). 

In order to have a unique solution we must require some additional conditions, 
known as entropy conditions. One of the most useful is the Lax condition, 
which says that a shock connecting the states u~ and , travelling with speed 
7 = Xi{u~,u~^) (Ai(u~,M+) is an eigenvalue of the averaged matrix 
cf. [3]) is admissible if 

K{u-)>\{u~,u+)>\{u+). (13) 

Now we define two types of curves: fix a state uq G M". Let ri{u) be the i-th. 
eigenvector of A{u). The i-th rarefaction curve through uq is the integral curve 
of the vector field and is denoted by cr i— > Ri{a){uQ). The i-th shock curve 
through Uq is the curve of states u which can be connected to the right of uq 
by an i-shock, satisfying the Rankine-Hugoniot conditions. It is denoted by 
a I— *■ Si{a){uo). The i-th rarefaction and shock curve are tangent to the ri(u) 
at Uq. 

Next we study the Riemann problem 



ut -\- f{u)x = 

x<0, (14) 



u(0, x) 



w+, X > 0. 



Under the assumption that the system is strictly hyperbolic with smooth coeffi- 
cients and each i-th characteristic field is either genuinely nonlinear or linearly 
degenerate, there exist three special cases: 

(1) Centered rarefaction waves: the i-th characteristic field is genuinely non- 
linear and w+ lies on the positive i-rarefaction curve through u~ , i.e. 

= Ri{a){u^) for some a > 0. Then the solution of is the centered 
rarefaction wave: 

{u" , X < tXi{u^), 

Ri{s){u~), x — tXi{s), s e [0, cr] 

U+, X > tAi(M+). 
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(2) Shocks: again the i-th. characteristic field is genuinely nonlinear, but this 
time the state u"*" lies on the i-th shock curve through u~ , i.e. u"*" = 
Si{a){u~). Denote the Rankine-Hugoniot speed of the shock Ai(u^,M+) 
by A. Then the solution of (HH) is the shock 



(3) Contact discontinuities: the «-th characteristic field is linearly degenerate 
and lies on the i-th rarefaction curve through u^, i.e. w+ = Ri{a)(u~) 
for some tr. Then the function is again a solution, but this time called 
the contact discontinuity. 

The parameter a, for which u+ — Ri{a){u^) or m+ = S'i(a)(u~), is called the 
wave strength. 

Therefore, if u+ lies on the rarefaction or shock curve the solution of the Rie- 
mann problem (|14p is one of the elementary waves - a centered rarefaction, a 
shock or a contact discontinuity. Otherwise, for m+ sufficiently close to u^, the 
Riemann problem (|14p can be decomposed in n auxiliary Riemann problems, 
which can be solved by an elementary wave. Piecing together those solutions 
we obtain a solution of the initial Riemann problem (jl4[) . 

2 Symmetry Groups of the Systems of Conser- 
vation Laws 

After we gave the brief overview of the notation and results from symmetry 
group analysis (classical and in the generalized setting) and conservation laws, 
we turn our attention to concrete systems of conservation laws. As we have 
just seen, a system of conservation laws is a system of first order partial dif- 
ferential equations. We introduced symmetry groups of system of differential 
equations as local transformation groups which act on the space of independent 
and dependent variables, transforming the solution of the system to other solu- 
tions. Also, we defined associated and fc-strongly associated symmetry groups. 
The aim of this section is to verify the results given in the introduction in two 
examples of conservation laws. 

2.1 A Model System of Two Strictly Hyperbolic Laws 

The first system we consider is 




X < Xt, 
X > At. 



(15) 



Ut + {u^ - v)x =0 




(16) 
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with initial conditions: 

»(-°)-{:; III "<-»)={::: 

A motivation for studying this system arises from some physical models like a 
model for a nonlinear elastic system or a model for the evolution of ion-acoustic 
waves. 

We start by calculating symmetry groups of the system (|16p . using the procedure 
described in the introduction. 

(1) Since v denotes one of the dependent variables, denote the infinitesimal 
generator by w: 

w = ^(x, t, u, v)dx + t{x, t, u, v)dt + 4>{x, t, u, v)du + "^l^ix, t, u, v)dy. 



(2) From ^ we get the first prolongation of this vector field 
pr(i)w = w + 4>''du^ + 4>'du, + Vdv^ + ^p'd, 
and by ^ we calculate 



= 0a; - 




\- (j)yVx 


^x'^x 


^u^x 


^v^x'^x 


-TxUt 


- TuUxUt 


- TyUtVx 




H 4'uUt ' 


t- (t)vVt 


- ^tUx 


- iuUxUt 


- ivUxVt 


~ TtUt 


- TuUl 


- TvUtVt 




1- ll^uUx - 


\-1pvVx 




^u'^x'^x 


- iwi 


- TxVt 


- TuUxVt 


- TyVxVt 


i/j* = Vt - 


\- IpuUt - 


\- ^vVt 


- £.tVx 


- ^uUtVx 


- ivVxVt 


- TtVt 


- TuUtVt 


- TvV^ 



(18) 

(3) Now we have 

Ai{x, t,U,V,Ux,Vx,Ut,Vt) = Ut + 2uux - Vx 
A2{x, t,U,V,Ux,Vx,Ut,Vt) =Vt+ U^Ux - Ux, 

therefore we need to solve the system 

pr(i)w(Ai) + 2(j)Ux + 2(lfu - V'^ = 
pr(i)w(A2) = + 2(j)uux + (jy^'u^ - = 0, 

whenever ut — ~2uUx + Vx and vt — —u^Ux + Ux- Inserting psp in these 
equations and replacing ut by —2uUx + Vx and Vt by —v?Ux + Ux, whenever they 
appear, we arrive at 

(t>t + 4>u{-2uUx + Vx) + 4>v{-u'^Ux + Ux) - itUx - ^uUx{-2uUx + Vx) 
-S.vUxi-U^Ux + Ux) - Tt{-2uUx + Vx) - Tu{-2uUx + Vx)"^ 
-Ty{~2uUx + Vx){-u'^Ux + Ux) + 2(t>Ux + 2u[(j)x + (j)uUx + (t>vVx 
"i.xUx - ^uul - ^vUxVx - Tx{-2uUx + Vx) ~ Tu,Ux{~2uUx + Vx) 
-Tyi-2uUx + Vx)Vx] - [ipx + IpuUx + IpvVx " £,xVx " £.iiUxVx 
-ivvl - Tx{-U^Ux + Ux) - TuUx{-U^Ux + Ux) - TyVx{-u'^Ux + Ux)] 

= 



9 



-CvVxi-y-'^Ux + Ux) - Tt^-U^Ux + Ux) - Tu{-2uUx + Vx){-U^Ux + Ux) 
-Tyi-U^Ux + UxY + 24)UUx + (u^ - l)[(t>x + 4>uUx + (t>vVx - £.xUx - ^uul 
-£.vUxVx - Tx{-2uUx + Vx) - TuUxi-2uUx + Vx) - Ty{-2uUx + Vx)Vx] 

= 0. 

(4) Apparently, solving this system is quite complicated. Hence, we are going 
to look only for projcctable symmetry groups. So, assume that ^ and r only 
depend on x and t. Then we have 

4>t - 2(j)uUUx + (t>uVx - 4>vu'^Ux4'vUx - £,tUx + 2TtUUx - nVx + 2(j)Ux 

+ 2(j)xU + 2(j)uUUx + 2(j)yUVx - 2^xUUx + ATxU^Ux — 2TxUVx -tpx - 'ipuUx 

- fpvVx + ^xVx - TxV?Ux + TxUx = 

tpt - 2lpuUUx + i'uVx - tpvU^Ux + IpvUx - £,tVx + TtU^Ux - TtUx + 24>UUx 
+ ipxU^ + (t>uV^Ux + (j^vU^Vx - S.xu'^Ux + 2TxU^Ux - TxU^Vx - (t>x - 4>uUx 

- (f^vVx + £.xUx - 2TxUUx + TxVx = 0. 

These equations are in fact polynomials of free variables x,t,u,v,Ux and Vx ■ The 
solution will be found by looking at their coefficients on the left and right hand 
side of the equations. Since the functions 0, ip and their derivatives depend on 
x,t,u and v we equate the coefficients of 1, Ux and to 0. Then we arrive at 
the following equivalent system: 

4>t-i'x + 2u(f>x = 

2(j) + (j)y - - + Tx + u{2Tt - 2^x) + U^i^Tx -(l)y) = 

K - i-'v + Cx -n + u(2(f)y - 2tx) = 

Ipt- 4>x+ u'^'t'x = 

i^v - 4>u - Tt + ^x + u{24) - 2%py, - 2Tx) + U^{(t>u - ■01) - Ca; + n) + 2u^Tx = 
■^u - <t>v - £,t + Tx + u'^{<j)v -Tt) =0 

(5) The general solution of this system is: 

^{x,t) = CiX + est + C4 
T{x,t) = Cit + C2 
(f){x,t,U,v) = C3 
1p{x,t,U,v) = C3U + C5 

ci — C5 arc arbitrary constants. The linearly independent infinitesimal generators 
of projectable symmetry groups are: 

wi = xdx + tdt 
W2 = dt 

W3 = tdx + du + udy 
W4 = dx 
W5 = dy 
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(6) It remains to compute the corresponding one-parameter symmetry groups. 
Hence, the one-parameter group Gi, generated by the vector field Wi, is a 
sohition to the system of ODEs 

x(r]) — x{r]) 

with initial data a;(0) = x and t(0) = t, i.e. 

grj : (x, t, u, v) (e^x, e^t, u, v) . 

Since Gi is a symmetry group, ^ implies that if u and v are solutions of (|16p 
so are the functions 

u : {x,t) ^ u (e^''x, e^^t^ 
V : {x,t) ^v{e-~'^x,e~'^t) . 

We repeat the same procedure for the remaining symmetry groups and calculate 
that G2, G4 and G5 are translations of i, x and v, respectively. Finally, the action 
of G3 is given by 

/ 772 \ 

g,f : {x,t,u,v) —> lx + r]t,t,u + ri,v + riu+—j, 

and the functions 

u : {x, t) u(x — rjt, t) + r/ 
V : {x, t) v{x — rjt, t) -|- riu{x — r]t, t) H — ^ 
are solutions of the system whenever u and v are. 

Therefore, we calculated all projectable symmetry groups of the system (|16p 
and all transformed solutions. As we saw, the calculation of non-projectable 
symmetry groups is rather complicated on the one hand, and on the other, as 
was mentioned in the introduction, it is enough to study projectable groups if 
the solution is in Q, V or if it is a solution in the sense of association. We 
recall from |B] that the system has a solution in the sense of association: 
the Jacobian matrix of (11611 is 

2u -1 
-1 



A = Df 



the eigenvalues are Ai(u, v) = u — 1 and \2{u,v) = u-l- 1, and the corresponding 
right eigenvectors are ri{u,v) = [1 u + ly and r2{u,v) = [1 u — 1]"^. Since 
DXi{u,v) ■ ri{u,v) > 0, i — 1,2, it follows that both characteristic fields are 
genuinely nonlinear and hence the solution consists only of centered rarefaction 
waves and shocks. 
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The rarefaction curves are calculated as the integral curves of the vector fields 
ri and 

Ri ~ {{u,v) : V — -u^ + u + ci} 

i?2 = {{u, v) -.v = ^u^ -u + C2}, 
and the shocks are found from the Rankine-Hugoniot equations: 



«-^;o^(^,-^,o)(^^T^l-^^4?^)' for|"-"o|<12. (19) 
The corresponding shock speeds are 



7 = "0 H — ± \/ 1 



2 V 12 ' 

where the sign — refers to the speed of 1-shock, and + to the speed of 2-shock. 
The Riemann problem p6|) - p?)) has a classical solution for each {u,v) lying in 
the area bounded by 

J{u,v) = u) : u = + 7i + ^ + Wo - - uo A u > Wo - sj 

Ji{u,v) — ^(u. v) : (u. v) satisfies A u < un — sj 

J2{u,v) = w) : w = - u - ^ + Wo - i^o + A u > Mo - sj 

The remaining {u, v) are in the exterior of this area, which we denote by Q and 
which is divided by the curves 

D{u, v) = {{u,v) : V = vq + + (1 — uq)u — uq /\ u < uq — 3} 
E{u, v) = v) : V ~ vq + {u — uq)(uo — 1) A u < uo — 3} 

into three open regions. In each of them the solution consists of a singular shock 
wave which is given by: 

U(x, t) = G{x — ct) + si(t){aQd^ (x — ct) + aid^{x — ct)) 

V{x, t) = H(x - ct) + +S2{t){l3fiD-{x - ct) + (3iD+{x - ct)), ^'^^^ 

where G{x — ct) and H{x~ct) are generalized step functions (cf. [13], Def. 1(a)), 
D{x — ct) — (3oD~{x — ct) + PiD~^{x — ct) is an S(5-function with value (/3o, /9i), 
Po+ Pi ^ I (cf. [T3I, Def. 1(b)), d{x - ct) = aod-{x - ct) + aid+{x - ct) is 
an 3'SD-function with value (ao,Q;i) (cf. [13], Def. 3 with Ex. (ii)), such that 
D(x — ct) and d{x — ct) are compatible and 

~c[G] + [G'] - [H] - 0, (21) 
S2{t)^slit){al + al), (22) 

■S2{t) = (Tit, (71 = c[H] - i[G3] + [G], di > 0, (23) 

cs2(<) = sl(t){aQUo + aiui). (24) 
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The function S2{t) is called the strength of the singular shock wave and is the 
most important part of the solution which has to be uniquely determined, ao 
and ai can be chosen such that + af = 1, hence the condition ([^ becomes 

S2{t) = slit), 

and the condition 

UqUo + a\ui = c. 

We are going to show that this is a solution in the sense of 1-strong association. 

2.1 Theorem. The solution i20\) of the system il6]) is a l-strongly associated 
solution to hlb]). 

Proof. In order to show this we use Definition [Lj] Let i? be a bounded subset 
of Cl{R X [0,00)). This means that there exists K CC K x [0,cx)), with the 
property supp Q K for each tp d B and with 

sup {Id^ifiix, t)\ : (p e B, \a\ <1} < CO . 

{x,t)eK 

It suffices to show that there exist representatives and of the solution (|20p 
such that 

({UMx,t) + (U^e~ye)x{x,t))<p{x,t)dxdt\=Q, (25) 



lim sup 



and 



lim sup 

^^Oijses Jrx[o,oo) 



({V,)t{x,t) + {^U^^-UMx,t)^ip{x,t)dxdt\=Q. (26) 
Look first at Let ip ^B. Then 

((;7,(x,t))^+(c/2(x,t)-K(a^,t)^ 



/ 



Rx[0,oo) 



(y9(a;, t) dx dt 



Ge{x — ct) + si{t){aod {x — ct) + aid'^ix — ct)) 

x[0,oo) 

+ |Gg(a; — ct) + 2si{t){aoUod~ [x — ct) + aiuid'^{x — ct)) 
+sl{t){al{d^)'^{x - ct) + al{d+)^{x - ct)) - He{x - ct) 



-S2{t){PoD-{x - ct) + PiD+{x - ct)) 



(p{x.t) dx dt 



cdxGe{x — ct) + s'i{t){ai)d [x — ct) + aid'^{x — ct)) 

x[0,oo) 

—csi(t)dx(aQd~ ix — ct) + aid'^{x — ct)) 
+dxGe{x — ct) + si{t)dx{ctaUod^ {x — ct) + aiuid'^{x — ct)) 
+sl{t)dx{al{d-)^{x - ct) + al{d+)^{x - ct)) 



-dxH^ix ~ ct) - S2{t)dx{l3()D- {x ~ ct) + (3iD+{x - ct)) 



ip{x, t) dx dt 
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Ix [0,00) 



s'i{t){aod {x — ct) + aid'^{x — ct)) ip{x, t) dx dt 



(1) 



— cG^{x ~ ct) — csi{t){aod {x — ct) + aid'^{x — ct)) 




+ G^(x — ct) + si(t){aQUod {x — ct) + aiuic? (x — ct)) 

(4) (5) 

+ sl{t){al{d-)^{x - ct) + al{d+)^{x ~ ct)) 



(6) 

He{x - ct) - S2{t){l3oD- {x - ct) + (3iD+{x - ct)) 
It) 



(Px{x, t) dx dt 



(8) 



(*) 



For the first member of tliis sum we fiave 

s'i{t){aod^ (x — ct) + aid'^{x — ct))Lp{x, t) dx dt 



Rx[0,oo 
oc 



+"1 



/I /x — ci — 4e\ 1 /x — ct — 6e\ \ 1/2 



-X + ct — 6e\ \ 1/2 



) - M- 



(p{x, t) dx dt 



Here we used tlie 3'SD-fimction from [T3]. Introducing suitable substitutions 
we obtain 

[ s[{t) / ao\/e y-(l>{z) j ip{—ez~'4e + ct,t)dzdt 
f°° f°° (\ \i/2 

' s'i(^) y '^oV£'(^2'^('^)j ^{—ez — ^e^ct,t)dzdt 



I dz dt 



f°° /I \i/2 
[{t) / ai^/ei^-(l){z)j ip{ez-4e + ct,t)i 

r° f°° (\ \i/2 
' y '^iv^(^2*^^^v ^{sz — &£ + ct,t) dz dt. 



Applying the Lebesgue dominated convergence theorem two times successively 
to the corresponding sequences we conclude that this term tends to as e ^ 0. 
A similar argument shows that each of the terms in the sum with the functions 
d^ or (d^)^, i.e. (3) and (5), also goes to as e ^ 0. So, look now at (2). We 
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have 



—cGe{x — ct)ipx{x, t) dxdt 

Tlx [0,oc) 

oo pct~e roc POO 

/ cuo(px{x,t) dx dt + / / cuiipx{x,t) dx dt 

J-oc JQ Jct+e 



cuq / ip{ct — e,t) dt — cui / ip{ct + e,t)dt 



(cuq — cui) I (p{ct, t) dt 



/•oo 

= -c[G] / ifi{ct,t)dt, 
Jo 



where we again applied the Lebesgue dominated convergence theorem. For (8) 
we obtain 



/ -S2{t){f3oD-{x - ct) + PiD+ix - ct))ifxix,t) dxdt 

S2\t) 01 fx(X, t) dx dt 

J-oo £ ^ e ^ 

°° -"^ , jx-ct-2e\ , , , , 
(t) — 0^ j(px[x, t) dx dt 

1 

S2{t) I f3Q4}{z)ipx{ez — 2s + ct,t) dz dt 



S2[ 

Jct+e 



+ / S2(^) j l3i(j){z)ipx{ez + 2e + ct,t)dzdt 

/>! poc 

S2[t)(3o(px{ct,t) I 4'{z)dzdt+ / S2{t)Pi(px{ct,t) I if){z)dzdt 



(/3o+/3i) / S2{t)^x{ct,t)dt 

Jo 

S2{t)ipx{ct, t) dt. 



Repeating this for the remaining terms yields that (4) tends to [G^] t/5(ct, t) dt, 
(6) to {al + al) sl{t)(px{ct,t) dt, and (7) to ~[H] (p{ct,t) dt, as e ^ 0. 
Therefore, 



(*) —i -c[G] / ip{ct, t) dt + [G^] / if{ct, t) dt 
Jo Jo 

/•oo poo 

+ {al + al) sl{t)(p{ct,t)dt - [H] Lp{ct,t)dt-I S2{t)<fi{ct,t) dt 
Jo Jo 

= 0, 



by (|2ip - ((24l) and ((25|) is satisfied. Similarly we conclude that it is also true for 
Hence, the solution (|20p is a 1-strongly associated solution to ^TE\i . □ 
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The projectable symmetry groups calculated at the beginning of this section 
transform 1-strongly associated solutions to to other 1-strongly associated 
solutions, as shown by the following 

2.2 Theorem. The symmetry groups Gi - of the system U6\) are -45^- 
symmetry groups. 



Proof. By Theorem 11.51 it suffices to show that Gi - Gs are slowly increasing 
and have a factorization (jlOp such that Q depends only on rj and x. First we 
consider Gi. The action of Gi is given by 

gjj : (x, t, u, v) — > {e^x, e^t, u, v) . 

Since is the identity it follows that the map 

{u,v) ^ ^g{x,t,U,v) 

is slowly increasing, uniformly for x and t in compact sets. It is easy to see that 
this is also true for the remaining groups. Next, for Gi we have 



Ai(e-''x, e-^t, pr(i)u(e-%, e^^i), pr^^^'y(e-"a;, e^'^i)) = er^ut + 2e-'^uu^ 



= e-"Ai 

— e^'^Vt + e^^u^Ux — e~' 
= e-^A,, 



where Ai and A, denote the first, respectively the second equation of the system 
(|16p . The matrix form of this factorization is given by 



" Ai ■ 




' e-'^ ■ 




■ Ai " 


. ^2 . 




e-'' 




. ^2 . 



Therefore, the matrix Q depends only on ry and Theorem 11.51 provides that the 
Gi is ^5^-symmetry group. The factorizations for G2-G5 are 



G2, G4, G, 



5 ■ 



" Ai " 




" 1 


" 




" Ai " 


. ^2 _ 







1 




. ^2 . 


" Ai ' 




" 1 


" 




" Ai " 


A2 






1 




. ^2 . 



G, 



hence again we conclude that these four symmetry groups are also ^5^-symmetry 
groups. □ 
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2.2 Zero Pressure Gas Dynamics Model 

The next system of conservation laws we consider is given by 



Ut + (uv)„ — 

(27) 

{uv)t + {uv'^)x = 

with the same initial conditions p7p as in the previous case. This Riemann 
problem is a zero pressure gas dynamics model, where u is a density, hence 
nonnegative, and u is a velocity. 

The quasilinear form of (P7|) is obtained by differentiating: 



Uf + UrrV + UV,r = 

(28) 

UtV + UVt + UxV + 2uVVx — 0. 



If we compute Ut from the first equation of (j28p and insert into the second one 
we arrive to the following system: 



Ut + UxV + UVx =0 , , 

X ^29) 

u{vt + Wx) = 0. 

From the second equation it can be seen that one possible solution is w = 0, i.e. 
vacuum state. Therefore, we consider the other possibility Vt + vvx = 0, looking 
at the system 



Ut + UxV + UVx = ^^^^ 
Vt + VVx = 0. 

The eigenvalues of this system are Ai(u, u) = A2(m, w) = w, thus the system is 
weakly hyperbolic. The corresponding right eigenvector is r{u,v) = [0 1]"^, so 
both characteristic fields are linearly degenerative. Hence, only contact discon- 
tinuities can appear as a solution and we calculate them: let x = ct he a curve 
of discontinuities of the system ([29|l . Along this curve the Rankine-Hugoniot 
conditions must hold: 

c[u] = [uv] 
c[uv] = [uv'^]. 

Equating c from these equations yields 

= [uv]'^, 

so 

uoUi{vi - vo)"^ = 0. 
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Therefore, there exist three solutions: uq — 0, ui — and vq = vi. For uq = 
we calculate c — vi, Xi{uo,vo) — and \i{u\^v\) = wi, i = 1,2 and similarly 
for u\ = 0. In the third case c = ~ v\ and also Ai(Mo,t'o) = = 
vq = vi, thus the Lax entropy condition is satisfied. Hence the initial 
conditions (uq^vq) and can be connected by contact discontinuity only 

when vq — vi. Finally, combining contact discontinuities and vacuum states we 
obtain the classical solution of the Riemann problem (j27p when vq < vi: 



{u,v){x,t) = { {0,v{x,t)), vo<j<vi 
{ui,vi), f > Vi. 

In the case when vq > vi this solution is not uniquely defined and certain 
nonregularities appear, which is studied in detail in [13] . In that case the solution 
of (|27p is again a singular shock wave 

U{x, t) = G{x ~ ct) + si{t)iaoD- + aiD+) + S2(t)(/3od" + Pid+) 

V{x, t) - H{x - ct) + S3(t)(7od" + lid+), ^^^^ 

where G and H are generalized step functions, D and d are compatible SJ- and 
3SD-functions (cf. [13 , Def. 3 with Ex. (i)) and 

si{t) = (7i<, (71 = c[G] - [GH], (71 > 0, (32) 

ao = , ai = , (33) 

aiiaovo + aiwi) (7ic c[GH] - [GH'^], (34) 
-S2{t)sUt) ^ si{t), (35) 
ao(wo - cvq) + ai{vl - cvi) = Polo + /'i7i- (36) 

This time the function si{t) denotes the strength of the singular shock wave. 
As for the first system (|16p we are going to show that this solution is also a 
1-strongly associated solution to (P7)l . 

2.3 Theorem. The solution 131]) of the system P7| ) is a 1-strongly associated 
solution to [2T 



Proof. The proof is similar to that for the system so we have to show 

that there exist representatives 11^ and of the solutions U and V defined in 
(|3T|) . such that for arbitrary set B C C;?°(M x [0, oo)) bounded in C^i(R x [0, cx))) 
the following holds: 

limsupl/ ({U,)t{x,t) + {UeV,)^{x,t))^{x,t)dxdt\^0, (37) 

and for the second equation: 

limsupl/ ({U,V,)t{x,t) + {U,V,^),{x,t))^{x,t)dxdt\=0. (38) 
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Let B be a bounded subset of Cl(R x [0,oo)). First we prove (P7|) . Let ip € B. 
Then 



Rx[0,< 



>f{x, t) dx dt 



Ge{x - ct) + si{t){aoD~{x - ct) + aiD+{x - ct)) 

Rx[0,! 

+S2it){Pod-{x - ct) + l3id+ {x - ct))\ + I Geix - ct)H,{x - ct) 



+si{t)s3{t){aoD {x-ct) + aiD+{x~ct)){-/ad^{x~ct)+^id+{x~ct)) 
+S2 (t) S3 (i) (/3o (cc - ct) + /3id+ (x - ct) ) (70 (x - ct) + 71 (i+ {x - ct) ) 
si{t){aoVQD^ {x — ct) + aiViD'^{x — ct)) 
+S2{t){f3oVod^ {x — ct) + /3iwid+(x — ct)) 



+S3it){'^oUod {x - ct) + 7i-uid+(a; - ct)) 



/ 



ip{x.t) dx dt 



cd^Ge{x - ct) + s[{t){aoD^{x - ct) + aiD"^{x - ct)) 



—csi{t)dx{anD {x — ct) + aiD^{x — ct)) 
+s'2{t){l3od- {x - ct) + I3id+{x ~ ct)) 
-cs2{t)dx{f3od~ {x - ct) + /3id+(a; - ct)) 
+dx{G,{x-ct)H,{x~ct)) 

+si{t)s-i{t)dx(^{aQD^ [x — ct) + aiD^{x — ct))(7od^(a; — ct) 

+7id+ (a; -ct))) 
+S2{t)s3{t)dx({Pnd-{x - ct)+Pid+{x - ct)){"fod-{x - ct) 

+71 d+ (a; -ct))) 

+si{t)dx{ctQVoD^ {x — ct) + aiViD^{x — ct)) 
+S2{t)dx{l3aVod^ {x - ct) + - ct)) 



+S3(t)dxi-^oUod {x - ct) + 7iuid+(x - ct)) 



ip{x, t) dx dt 



Rx[0,! 



s[{t){aoD~{x - ct) + aiD+{x - ct)) 



(1) 



s'2{t){Pod-{x - ct) + - ct)) 



(p{x, t) dx dt 



(2) 
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Rx[0,oo) 



cG^{x — ct) — csi{t){aoD {x — ct) + aiD^{x — ct)) 



(3) (4) 

- cs2{t){l3Qd~ {x - ct) + I3id^{x -ct)) + Ge{x - ct)H^{x - ct) 

v ' V ' 

(5) (6) 

+ si{t)s3{t){aoD~ {x-ct) + aiD+ {x-ct)){jod~ {x-ct)+-fid+ (x-ct)) 

(7) 

+ S2it)s3{t){Pod^{x-ct)+Pid+{x-ct)){-/od-{x-ct)+-fid+{x-ct)) 

(8) 

+ si{t)(aQVoD^ {x — ct) + aiViD^ {x — ct)) 





(9) 




+ S2{t)ij3aVod {x 


- ct) 4 


- (3iVid'^{x — ct)) 




(10) 




-\- S3{t){'JoUod (x 


-ct)4 


- ^\U\d^{x — ct)) 




(11) 





ipx{x, t) dx dt 



(*) 



Consider now each of the terms in the last sum, hke in the proof of Theorem 
2.11 For (1) we have: 

s'i{t){aoD^ {x — ct) + aiD^{x — ct))(p{x, t) dx dt 

Rx [0,oo) 

,,,/ao,/a; — ci + 2e\ ai ^ / x ~ ct — 2e\\ , ,, , 
iW(t'^( 1 j'^e'^K i )y{x,t)dxdt 

Next, spht up this integral as a sum of two integrals and take suitable substi- 
tutions. This yields 

ni 
s\ (t)aQ(^{z)ifi(ez — 2e + ct, t) dz dt 

noo 
s\{t)ai(t>{z)Lp{£z + 2e + c<, t) dz dt 

Here we apply the Lebesgue dominated convergence theorem, first to the se- 
quences ( 'p{£z — 2e + ct, t) dz)^ and ( (p{ez + 2e + ct, t) dz)^, and then 
to the sequences {(p{ez — 2e + ct,t))^ and {fiez + 2e + ct,t))^. Then the last 
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integral 



»1 

s'i{t)ao'f{ct,t) / (f>{z) dz dt + / s[{t)ai^p{ct,t) / <j){z)dzdt 



Jo J-oo Jo 

/"OO 

(ao + ai) / s'i(t)ip{ct,t) dt 







s'i{t)ip(ct, t) dt, 

since by definition of S^-functions, / = 1 on the domain of (j) (and 

that is the interval [—1, 1]) and ao + ai = 1 by assumption. It is obvious that a 



procedure similar to this one and those from the proof of Theorem l2.1l is repeated 
for each of the remaining 9 terms (some of them are explicitly calculated in the 
proof of Theorem 1 2 . 1 p . The only difference is that for the solution of the system 
(PT|) we used SSD-functions instead of 3'SD-functions from the solution of (fTB|). 
This implies that all terms with or (d^)^ tend to as e ^ 0. 
Joining all together we have 

/•OO POO /"OO 

(*) ^ / s'^{t)Lp{ct,t) dt ~ c[G] / ip{ct,t)dt+ I csi{t)ip^{ct,t)dt 

dt 



"'0 Jo 

POO POO 

+ [GH] / Lp{ct,t)dt~ I {uqVq + aiVi)si{t)Lpx{ct,t) 
Jo Jo 



= 0, 

since from (15^ si(^) — "^[G] — [GH] and aouo + aiVi — c from and the 
condition ag +ai — 1. Thus ((3T|) is a 1-strongly associated solution to the first 
equation of (P7|) . Analogously, it can be seen that the same holds also in the 
case of the second equation, which proves the claim. □ 



The next task is to calculate symmetry groups of ((27|) . Let us recall that (jST]) . 
((28|) and ([29]) are equivalent systems for all smooth solutions. Also, these sys- 
tems are equivalent in the Colombeau algebra, since the elements of this algebra 
are equivalence classes of sequences of smooth functions. Therefore we look for 
the symmetry groups of the quasilinear system (I29|) . We start with the following 

2.4 Theorem. Let A(a;, m^")) = and Ai(a;, u*^"^), i = \, . . . ,k, be nondegen- 
erate differential equations on M d X x U such that A can be written as a 
product 

k 

4=1 

If we denote the corresponding algebras of infinitesimal generators of symmetries 
of A and A^ by q and Qi respectively, i = 1, . . . , k, then 

k 

n 0^ ^ 0- (40) 

1=1 
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Proof. Let V G nJL^fli- Then v G g,, Vz = 1, . . . , k, i.e. v is a generator of a 
local one-parameter symmetry group of each equation 

A,(x, = 0, i = l,...,k. 

By the infinitesimal criterion and ^ we may write 

pr(")v(A,(x, =Q,-A,(x,7/(")), i = l,...J, (41) 

with well-defined functions Qi, i — 1, . . . ,1. Since pr'^"W is a vector field on the 
n-jet space M^"^ the Leibniz rule for the product derivative yields 

fc k 

pr(")v(A) = pr(")v( [] A,) = ^ Ai • . . . • pr(")v(A,) • . . . • A^ 

i=l z=l 

fe 

= ^ Ai • . . . • Q,A, • . . . • Afc = Q • A, 

i=l 

where Q = Qi + . . . + Qi. Another application of Theorem 11.21 provides that 
V G which proves the claim. □ 

According to this theorem, the intersection of the symmetry groups of the system 
(|30p with the symmetry groups of u = will provide symmetry groups of the 
system 1^^. (It should be noticed that due to the inclusion in (HD|) not all 
symmetry groups of pS)) will be obtained. However, Theorem 12.41 is of great 
help, since a direct computation of the symmetry groups of (|29p is a very difficult 
task.) The symmetry groups of w = can easily be calculated. Namely, by the 
infinitesimal criterion ([5]) it follows that the infinitesimal generators of u = 
are obtained as solutions of v(u) = </) = 0, whenever u = (we assumed here 
that v — £,{x, t, u, v)dx + t(x, t, u, v)dt -I- 4>(x, t, u, v)du + i^ix, t, u, v)dv). 
Now we follow the procedure for calculating symmetry groups of the system 
dSQl). 
(1) 

w — (^{x, t, u, v)dx + t{x, t, u, v)dt + 4>{x, t, u, v)du + "^Jix, t, u, v)dv. 

(2) The first prolongation is given by 

pr(i)w = w + 0"a„^ + -I- ^^a,, -I- 

where 0^ , 0* , and V'* are the same as in ^TE\\ . 

(3) Since the equations of this systems are 

Al(x,t, U,V,Ux,Vx., Ut, Vt) = Mf + UxV + UVx 
A2{x, t,U,V,Ux,Vx,Ut,Vt) = Vt +VVx, 

we have to solve 

pr(i)w(Ai) = (j/ + vcj)^ + <j)Ux + 4>Vx + — 
pr(i)w(A2) = -0* + "0^1 + vi'^ — 0. 
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Again we look only for the projectable symmetry groups. Inserting having 
in mind that the partial derivatives of ^ and r with respect to u and v vanish, 
and then substituting ut by —u^v — uvx and vt by —vvx, we obtain 

- 0„UMa; - 0«WWa: " 4>vVVx - ^tUx + nVUx + nUVx + 0a;W 
+ + 4>vVVx - £,xVUx + TxV^Ux + TxUVVx + ijjUx + (jjVx + tpxU 

+ ll^uUUx + ikvUVx - ixUVx - TxUVVx = 

■0t - ll^uVUx - IpuUVx - IpvVVx - £,tVx + nWx + ll^Vx + ll^xV 
+ ll^uVUx + V'l^W^^x - ^a: + TxV^Vx = 0. 

(4) Coefficients of f , Ux and u^; equating with yield the following equations 

4>t + Vlpx + U^x = 

-£.t + TtV - ixV + Txv"^ + + liV' = 

— (f>uU + TtU + TxUV + + IpvU — ^xU + TxUV = 

tpt + tpxV = 

= 

-IpuU - S.t + TtV + l/j - £,xV + TxUV = 0. 

(5) The solution is 

^{x, t) = Cl+ C^t + (C2 + cst)x + CsX^ 
t{x, t) ~ Cq + (C7 + Cst)t + (C4 + C5t)x 

4>{x, t, u, v) = ua{x, t, v) 

^P{X, t, U, v) = C3+C2V + 2c5XV + Cs{tv + x) — v{c7 + 2cst + C5x) — v'^{c4+C5t), 

where a is a function which depends on x, t and v and satisfies the equation 

cs + C5V + at+ vax _ p ,^2) 

v 

The eight constants ci — cg generate eight linearly independent infinitesimal 
generators of one-parameter projectable symmetry groups, while a{x,t,v) gen- 
erates an infinite-dimensional group. From (|42p we see that a{x,t,v) must de- 
pend on constants C5 and cg. It is also clear that the function a is not uniquely 
determined. Hence, in order to calculate the infinitesimal generators of the 
projectable symmetry groups we choose one possibility for a: 

a{x, t, v) = —C5X — cg,t + /3(f). 

Now we can write all infinitesimal generators: 

wi = dx 

W2 = xdx + vdy 
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W3 = tdx + 

W4 = xdt — v^dv 

W5 = x^dx + xtdt — xudu + [xv — tv'^)dv 

We = dt 

W7 = tdt — vdv 

wg = xtdx + t^dt — tudu + {x — tv)d^v 
w/3 = uf3{v)du. 

(6) The one-parameter transformation groups generated by the vector fields 



{x + 77, t, u, v) 
{e^x, t, u, e^v) 
{x + 7]t, t,u,v + rj) 

V 

x,t + rjx, u, 



Wi — Wg 


and W/3 are: 




: {x,t,u,v) 


G2 


: {x,t,u,v) 


G3 


: {x,t,u,v) 


G4 


: {x,t,u,v) 


G5 


: {x,t,u,v) 


Ge 


: {x,t,u,v) 


G7 


: {x,t,u,v) 


Gs 


: {x,t,u,v) 


Gp 


: {x,t,u,v) 



1 



t 



1 — rjx' \ — rjx 
(a;, t + rj,u, v) 
(a;, e^t, u, e~^v) 
X t 



riv 



, (1 - 'qx)u, 



1 — ri{x — tv) 



(43) 



1 



, (1 — rit)u, rjx + (\ — rjt)v 



r]t' 1 — Tyt ' 

Since each of the groups in (|43p is a symmetry group of the system pop. from 
([l} it follows that if u and v are solutions so are the functions 



(1) u : (a;, t) u{x — r/,t) v : {x, t) v{x — rj, t) 

(2) u: {x,t) u{e^'^x,t) v : {x,t) ^ e''v{e'''x,t) 

(3) u : (a;, t) u{x — rjt, t) v : {x, t) v{x — rji, t) + r/ 



(4) u : {x, t) — > u{x^ t — "qx) v : {x, t) 



v{x, t — r]x) 



1 + r]v{x, t — rjx) 
X t \ N / a; t 



u , (1 + r/a;)w 



1 + rjtv ' 



1 + rjx 1 + 77a; 

(6) u : (x, t) u(x, t — 7]) V : {x,t) ^ v{x, t — v) 

(7) u: lx,t) ^ u{x,e^'^t) v : {x,t) ^ e~'^v{x, e^'H) 

X t \ f X t 

, rix + V \ , 

^ ' ^ ' ^ 1+T]t ^ ^ l + Tjt 

{(3) u:{x,t)^e''^^'"\{x,t) v : {x,t) ^ v{x,t) 

From the remark given after Theorem l2.4l it follows that all calculated symmetry 
groups are also symmetry groups of the system (j29p : for infinitesimal generators 
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W5, Wg and W/j the coefficients of du are yq^w, and e^^^'"^u respectively, 

hence they vanish when u = 0, while for the rest = 0. 

The matrix factorizations of ([27]) with respect to symmetry groups in p3|) are: 



Gi, Gg 



G 



2 • 



Gs: 







■ 1 


" 




Ai " 


. ^2 . 







1 




A2 _ 


" Ai " 




" 1 







" Ai " 


. ^2 _ 







e'' 






" Ai " 




" 1 


■ 




■ Ai " 


A2 






1 




A2 



G4 



Ai 

A2 



l+2riv 



l — rfv 



Ai 

A2 



Gs 



Ai 

A2 



l + 27)fU 

(l+r;a;)^(l+jjt'u) 



l — r}tv 



(l+rix)(l+ritvY (l+r]x){l+ritv)^ 



Ai 

A2 



Gt 



" Ai ■ 




■ e-'' 


. ^2 . 




e-2') 



Ai 

A2 



Gs 



Ai 

A2 



(!+>)*)■■' 

'qx 





1 



Ai 

A2 



G, 



' Ai " 




. ^2 . 





(1 - 77^/3' (z;))e'''^('') r;^'(w)e'"^('') 
—rjuv 





■ Ai " 




. ^2 . 



Therefore, the matrix of factorization Q depends only on x, t and rj for all groups 
except for G4, G5 and Gp. In these three cases the factor Q depends also on v. 
Beside that, from the transformed solutions we see that the groups G1-G3 and 
Gq-Gs are slowly increasing, uniformly for x and t in compact sets. Thereby we 
have proved the next 

2.5 Theorem. The symmetry groups Gi, G2, G3, Gg, G7 and Gg of the sys- 
tem \2'T^ are 1-strongly associated symmetry groups, i.e. transform l-strongly 
associated solutions to |i?7|j to other 1-strongly associated solutions. 



The remaining three groups G4, G5 and G/3 are not ^5^-groups for two reasons: 
first the condition that the map 

{U,V) 1-^ ^g{x,t,U,v) 

is slowly increasing, uniformly for x and t in compact sets, does not hold globally. 
Second, the solution (I31|) does not belong to the algebra Qoo, which is necessary 
by Theorem 11.51 
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Under certain assumptions on the solution {u,v) defined in (PT|) . and also on 
the parameter 77, these problems can be avoided. Namely, if we assume that 
the function v is nonnegative and ry > (then instead of a group we consider 
a semigroup) the symmetry groups G4 and G5 become slowly increasing, while 
for Gg it should be supposed that P{v) is a function of L°°-log-type. The second 
condition from Theorem 1 1.51 (i) would be fulfilled if we assume that the solution 
{u,v) belongs to the algebra Goo- 
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